Effect of Intrinsic Curvature on Semiflexible Polymers 
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Recently many important biopolymers have been found to possess intrinsic curvature. Tubulin 
protofilaments in animal cells, FtsZ filaments in bacteria and double stranded DNA are examples. 
We examine how intrinsic curvature influence the conformational statistics of such polymers. We 
give exact results for the tangent-tangent spatial correlation function C(r) = (t(s) .t(s + r)) , both in 
two and three dimensions. Contrary to expectation, C(r) does not show any oscillatory behaviour, 
rather decays exponentially and the effective persistence length has strong length dependence for 
short polymers. We also compute the distribution function -P(R) of the end to end distance R and 
show how curved chains can be distinguihed from WLC using loop formation probability. 

PACS numbers: PACS : 87.15.He, 87.16.Ka, 36.20.Ey, 82.35.Pq 



A cell hosts variety of polymeric filaments. Tubulin 
filaments in eukaryotic organisms and FtsZ filaments in 
prokaryotic organisms, are important members of the cel- 
lular cytoskeleton. They also play important role during 
cell division. Therefore in vitro properties of these fil- 
aments, both in their isolated filament form as well as 
in bundled form, have been of great interest [H 0, 01]. 
One of the important mechanical properties which tubu- 
lin and FtsZ filaments share is their intrinsic curvature. 
The mechanisms by which microtubulesgenerate pulling 
force during chromosome segregation [l| and FtsZ-ring 
generates contractile radial force during cell constriction 
[J], exploit their intrinsic curvature. 

Interestingly biopolymers with intrinsic curvature and 
torsion have been known for long. Prominent exam- 
ples are the alpha helices in proteins. They have also 
been extensively modeled [5|. But polymers with only 
intrinsic curvature (but no intrinsic torsion), have not 
received much attention. This is perhaps because the he- 
lical structures are visually easy to detect compared to 
intrinsic curvature without torsion, except in two dimen- 
sions (2D). Intrinsic curvature is evident in electron mi- 
croscope Q and and atomic force microscope (AFM)Q 
pictures of FtsZ and AFM pictures DNA [||. As our 
analysis will reveal, in three dimensions (3D), unless both 
L/lp and l p /Ro are of order one, where Rq is the intrinsic 
radius of curvature of the polymer, L its contour length 
and l p its persistence length, it is difficult to differentiate 
them, visually, from a WLC. In this letter we will focus 
on conformational properties of curved chains, that may 
distinguish it from a WLC. Such diagnostic markers are 
important in the context of ever increasing number of 
new biopolymers, with properties which cannot be ex- 
plained by the standard models like Rouse, Zimm, freely 
jointed chain (FJC)[3 or WLC, especially for short coun- 
tour lengths. We ignore excluded volume effects and self 
crossings which are negligible if L/ Rq < 1; but important 
otherwise, as well as for non-dilute polymer solutions. In- 
terplay of the additional length scale Rq with l p and L is 
of interest here. 

Equilibrium, conformational statistics of such curved 



polymers have been studied using Monte-Carlo simula- 
tion Q and auxiliary field theory Their dynamics 
have been studied by numerical solution of their nonlin- 
ear equation of motion 10]. The field theoretic calcula- 
tion [9( assumed zero average torsion for the curved poly- 
mer and used a spatially varying auxiliary field. They 
predicted an oscillatory decay C(r) ~ e~ r / lp cos(r/i?o) 
for the correlation function in three dimensions. Here 
lp is the effective persistence length. This oscillatory de- 
cay is geometrically suggestive, as the polymer may wind 
around itself (loop) due to its intrinsic curvature. This 
also implies that at large l p , for L > 2ttRo the peak of the 
distribution function P(R) of the end to end distance R 
will alternate between R = and R = 2Ro, correspond- 
ing to complete loops and half loops, respectively, as the 
length L is increased. 

We compute C(r) exactly both in 2D and 3D and 
show that the decay is purely exponential, while oscil- 
lations can be recovered in 2D through extra constraint 
on the sign of the preferred inter-bond angle and through 
imposing torsion in 3D [5j. The absence of oscillations 
can be understood as follows. Although intrinsic cur- 
vature constrains the magnitude of the local curvature 

~ Rq 1 , the direction of the curvature vector 
remain uncorrelated along the contour, leading to expo- 
nential decorrelation. In other words the local plane of 
curvature keeps changing randomly which disfavors for- 
mation of planar or helical loops. Although the func- 
tional form of C(r) cannot distinguish WLC from curved 
chains, the effective persistence length l p turns out to be 
N dependent. We also compute the distribution func- 
tion -P(R) for curved chains which shows difference with 
WLC. 

A polymer of length L is described by a space 
curve R (s) , s G [0,£]. The Hamiltonian of an in- 
trinsically curved polymer is given by [10], k^t ~ 

~2 Js=o [ K ( s ) ~ Rq 1 ] 2 ds. The curvature vector k(s) — 
^[^■, where i(s) is the local tangent to the curve (i(s) — 
anc l r is the intrinsic radius of curvature. Chang- 
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ing variable x = — , KgT — , 

This reveals that l p /L and I p /Rq are the two important 
dimensionless ratios of the problem. In the discrete limit 
the monomers are given by the position vectors Ri and 
the normalised bond vectors are U = (R^+i — Rj)/6, 
where b is the bond length. The discretized Hamiltonian 
reads 
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where A = Since the hamiltonian depends only 

on nearest neighbor bond angles, we can use a stan- 



dard property [HI, {U.i k ) = (i i .tVi) |fc_i| ■ This yields 
(t(s).t(s + r)) = cxp(— r/lp), where the spatial sepa- 
ration r — b\k — i\. Thus quite counterintutively, in 
3D, C(r) does not show any oscillation, but decays as 
cxp(— r/lp). We identify the effective persistence length 

. « « « f dn z cos8 ie - h i( e i) 

as l p = -b/ln(ti.U+i) and {U-U+x) = j dn . e - hM) — • 

Here 9i is the angle between hi and tj+i and the integra- 
tion is done by holding ti as the polar axis and integrat- 
ing ti+\ over the solid angle . The thermal weight is 

eM-hi), where h^) = [ V1 ^ - ^g] 2 (see EqU). 

(U.i i+ i) is given by 
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where y = V / Ll p /2N/R , yi = y / 2Nl p /L. Using this 
expression we plot l p /l p in FigfTJas a function of i?o an d 
N. At finite Ro,L, but large AT, 2 P /Z P — > 1. But inter- 
estingly the discreteness of the chain has practical im- 
portance. For example, DNA has a finite bond length 
approximately equal to the size of a base pair (bp) and 
experiments with short DNA strands (15 — 90bp) have re- 
vealed fl2f an apparent persistence length which is 3 — 4 
times lesser than the standard value ~ 50nm (150bp). To 
explain this Ref 12] has ivoked new physics at the scale 
of base pairs. Models with long range correlation in the 
intrinsic curvature disorder in base pairs has also been 
proposed [||. But our calculation shows that a simpler 
factor, such as an uniform curvature can lead to substan- 
tial decrease in l p . FigQ] shows that for low Rq and high 
l p the reduction in l p could be two folds. 

Since most curved polymers have been detected on 
2D substrates, 2D requires special attention. We will 
now compute C(r) in 2D. Here the angle between the 
bond vectors U and ffc, 6i t k can be written as a sum 
of the intermediate angles between the successive bond 
vectors. Denoting fc = i + I, 9i, i+ i = YljZi Note 



that the sign of the angles are important here and 

is defined as the angle the tj + i vector has to rotate with 

respect to the tj vector, and thus C(r) is given by 
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Note that if the symmetry of h(9) with respect to ±9 
is not broken then the integrals J_ are real and hence 
(cos0j,j+z) will decay exponentially, as in 3D. The inte- 
gral can be estimated by noting that the numerator gets 
its maximum contribution from 9 — 9 (and not — 9 , 
to break the symmetry) where the exponent [ ^ 1 ^^ 0Q - — 

is zero. In the continuum limit (b/R Q <C 1) the 
maximum occurs at a small value of 9q ~ b/Ro which 
allows the approximation cos 9 ~ 1 — . The result- 
ing integral can be evaluated analytically (Eqf4]) and its 
agreement with the numerical evaluation of the exact in- 
tegral in Eq|3]is excellent (figure not shown). 
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where r — l(L/N). For finite Rq,L, but large N, the A-dependent factor Re[..] — ► 1 and we get C(r) 




FIG. 1: The effective persistence length l p /l p with respect to 
R and N, at l p /L = 0.2 (left) and l p /L = 10 (right). 
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e~ r / l p cos(r/R ). Thus because of the preferred direc- 
tion of the angle 9 (as opposed to ±#o) the decay is 
oscillatory. This is analogous to imposing an intrinsic 
torsion in 3D. 
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FIG. 2: P(P) vs R/L for stiff chains (L/lp = 1) for different 
radii of curvatures, A) Rq = 0.1, B) Ro = 0.3 and C) Ro = 10 
(in units of l p ). At large Ro/l p the WLC limit (symbols) is 
reached while at small Rq the peak of P(R) can go past the 
flexible limit to smaller values of R, and thus can attain a 
globule size smaller than that of flexible Freely jointed chain 
(FJC). We used the same N for all the chains. Inset: Typical 
configurations of a) FJC, b) WLC and c) intrinsically curved 
chain, (b) and (c), both have l p /L = 1 and (c) in addition 
has Ro — 0.2. These configurations were generated using MC 
simulation. 



FIG. 3: Probability of forming closed loops P(R = 0) ver- 
sus contour length. From top to bottom R = 1.5,3,7.5 and 
20, at fixed l p = 10. At small l p /R a (bottom-most) WLC 
is recovered. We fitted all the curves (solid lines) with the 
approximate formula derived by Shimada and Yamakawa 
P(R = 0) = (896.32/Z) 5 exp[-14.054/Z + 0.246Z], for WLC 
in the range I = L/l p < 10. Clearly for large l p /R (the top 
two plots) the fit fails. That provides us a diagnostic tool for 
distinguishing a curved polymer from a WLC. 
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normalization constant for j dRP (R) = 1. Following 
Ref . [l3l] . P(R) can be connected to the partition func- 
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tion Z(f) — J dti ■ ■ ■ J dtjq exp 

where / is an external field (analogous to external 
magnetic field in Ising model). This connection is 
made via the reduced probability distribution p (z) = 
J dR P (R) 6(Rs~z), which is the probability that 
one end of the polymer is fixed at R = 0, and the 
other end lies in a given z plane. Using the definition 
of -P(R) it turns out that the Laplace transform of p(z), 
PU)=5- L dz exp(/z/g p(z)=Z(f)/Z(f = 0). For 
computation purpose we converted this Laplace trans- 
form to a Fourier transform by choosing / purely imag- 
inary. Substituting the Hamiltonian in Z (/) we get 

matrices \]A 



■ Ei N ^P { IXi 1 ^ + ^ (*i + + 
} . Z (/) can be computed using transfer 



We now compute the distribution function P(R), us- 
ing transfer matrix methods. We also verified these re- 
sults using a Monte-Carlo (MC) simulation which gave 
additional insights into the typical configurations of 
the curved chain (see inset of Figf2]). The distribu- 
tion function of the end to end vector R is P (R) = 



where the elements of the transfer matrix V are 
{ii\V\i i+ i) = exp[-hi + ^(U + U+i) • z]. After com- 
puting Z(f)/Z(f = 0) numerically, we inverse Fourier 
transform it to obtain p{z) and using the relation P(R) — 

~2nz ^dz" 1 \ z=R ( wmcu is obtained in Ref. 13] using to- 
mographic method and the isotropy of P(R)) we obtain 
P(R). 
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For semifiexible polymers with intrinsic curvature the 
energy of a configuration is proportional to l p , at fixed L. 
Therefore at large l p , little deviation of local R from Rq 
makes the energy cost for such fluctuations large. Such 
rare but important small angle fluctuations, which re- 
quires large sampling in a Monte-carlo simulation, can 
be efficiently captured by transfer matrix method. In 
this stiff limit the peak of the end-to-end distribution 
function P(R) is pushed away from the entropy domi- 
nated Gaussian peak (see Fig(2]) and forces the polymer 
to make circular arcs of radius Rq. But at small Ro, very 
small globules, smaller than even that of an FJC (of same 
N), can form (Figj2]) . Whereas for small value of l p the 
behaviour will be dominated by the flexible limit and the 
peak of P(R) moves towards the gaussian peak. 

In Fig[3] we show how loop formation probability 
P(R = 0) changes with contour length, in different 
regimes of l p /R. Ref||[ had reported similar result using 
a Monte-carlo sampling. They found that as Rq increases 
the peak of P(R = 0) unexpectedly shifts to L values 
lesser than 2ttRq. This was argued as thermal softening 
of the polymer. But it is unclear why thermal fluctua- 
tions should favour larger bending over smaller bending 
(relative to curvature Rq 1 )- We argue that actually the 
shift occurs due to the competition between two different 
peaks in P(R = 0), and it also gives us a diagnostic tool 
to differentiate a curved chain from a WLC. At small 
l p /Ro i.e., in the WLC limit P(R = 0) has a peak at 
approximately L/l p ~ 3 which emerge from a competi- 
tion of bending energy and entropy Q . Whereas for large 
lp/Ra the curved chain has a peak at L = 2ttRq which is 
driven by energetics of intrinsic curvature. Notice that 
for the upper most plot (in FigJ3J) the peak is located at 
L/2ttRo = 1, while in the lower most plot it is located at 
L/lp ~ 3. We note that RefQ had fitted their data for 
P(R = 0) with a formula in Ref[f|, which however was 
derived for a curved polymer with nonzero torsion. 

Our transfer matrix calculation is different from the 
usual one in that we can not use periodic boundary condi- 
tion which allows one to approximate the partition func- 
tion as X^ax, where X max is the largest eigen value of 
the transfer matrix V. This is because we work in the 
regime l p ~ L where the effect of the boundary con- 
ditions can propagate deep inside the chain. In other 
words, our system size is not larger than the correla- 
tion length and hence strong finite size effects are ex- 
pected. In Eq|5J Z(f) is obtained as a weighted sum 
over all the matrix elements of V 1 ^" 1 . In order to eval- 
uate y*' 1 numerically, the angular space f2 = (#,</>) 
for a bond vector t{ is divided into n = nx x ni bins. 
Thus the matrix V has dimension n x n. Note that 
while obtaining V 1 ^^ 1 in Eq[5] a matrix multiplication 
V 2 = j dtV\t){i\V is numerically implemented through 
a discrete sum: [V 2 ]ik = J dQjV(Sli,£lj)V(Qj ,flk) — 
c^2j =1 Vij smOjdjiVik — c[VSV]ik- Here the matrix 



Sji = Sji sin 6j , with no summation over j intended, 
[151 ] arose from d£lj — sin9jddjd(f>j and c = d6d(j) = 
(w/ni) x (27r/n 2 ). 

Finally, in order to get V 1 ^^ 1 we needed a multiplica- 
tion of the type VSVSVS . . . V. Since for high l p and 
Rq, the main contribution to Z(f) comes from configu- 
rations where the inter-bond angles are small (and even 
smaller if N is increased), the angular discretization has 
to be dense enough to pick up the contribution from the 
small intcrbond angles. This is computationally expen- 
sive. So we chose length discretizations as N — 2 = 2 m so 
that using 0(log 2 iV) number of matrix multiplications 
we could obtain the matrix (VS) N ~ 2 . 

In summary, we showed by exact calculation that both 
in 2D and 3D, intrinsically curved polymers give expo- 
nentially decaying tangent correlation, C(r), as semi- 
flexible polymers. But the apparent persistence length 
l p is substantially lesser than l p for short chains, when 
the ratio I p /Rq is large. Contrary to physical expecta- 
tion C(r) does not have oscillatory decay with period 
L — 2ttRq, unless we impose torsion in 3D or fix the sign 
of the preferred interbond bending angle in 2D. Curved 
biofilaments of FtsZ or tubulin can be distinguished from 
WLCs' through loop formation probabilities, which has 
not been measured yet for curved polymers. 
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